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Abstract 

We present measurements of nonproportionality in the scintillation light yield of bismuth germanate 
(BGO) for gamma-rays with energies between 6 keV and 662 keV. The scintillation light was read out 
by avalanche photodiodes (APDs) with both the BGO crystals and APDs operated at a temperature 
of ~ 90 K. Data were obtained using radioisotope sources to illuminate both a single BGO crystal in a 
small test cryostat and a 12-element detector in a neutron radiative beta-decay experiment. In addition 
one datum was obtained in a 4.6 T magnetic field based on the bismuth K x-ray escape peak produced 
by a continuum of background gamma rays in this apparatus. These measurements and comparison to 
prior results were motivated by an experiment to study the radiative decay mode of the free neutron. 
The combination of data taken under different conditions yields a reasonably consistent picture for BGO 
nonproportionality that should be useful for researchers employing BGO detectors at low gamma ray 
energies. 
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1. Introduction 

For an ideal scintillator, the number of optical 
photons produced by a gamma ray is linearly pro¬ 
portional to the energy deposited in the scintilla¬ 
tor by the gamma ray. However, many scintillators 
exhibit nonproportionality in their response, espe¬ 
cially for gamma-ray energies below 100 keV [J. 
Nonproportionality is typically defined by the ra¬ 
tio of the scintillation yield for a particular energy 
deposit to the scintillation yield due to complete 
energy deposit for a 661.7 keV gamma ray from a 
137 Cs source. Measurements of nonproportional¬ 
ity in the response of bismuth germanate (BGO) 
in the following energy ranges have been reported: 
4.5 keV to 662 keV 0, 6 keV to 1332 keV [3], 14 
keV to 1330 keV [4], 60 keV to 662 keV [5], and 11 
keV to 100 keV pQ. We report new measurements 
in the energy range from 6 keV to 662 keV. 

Our measurements and comparison to prior re¬ 
sults were motivated by an experiment to study 
the radiative decay mode of the free neutron. We 
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reported the first observation of this process in a 
pilot experiment in which gamma-rays were de¬ 
tected by a single BGO crystal [6, 7 . To accu¬ 
rately test theoretical predictions for the photon 
energy spectrum, a higher precision experiment 
has been performed with a 12-element BGO de¬ 
tector [8 . Analysis of this experiment requires 
both Monte Carlo modeling of the energy deposit 
spectrum in the BGO crystals and understanding 
of the light yield from BGO. For our minimum 
detectable gamma ray energy deposit of ^5 keV, 
the nonproportionality for BGO is an important 
source of uncertainty. Hence we performed mea¬ 
surements to determine the nonproportionality ob¬ 
served for our crystals, temperature, electronics, 
and other operating conditions. Our results were 
determined by both tests on a single BGO crys¬ 
tal with radioactive sources in a small test cryo¬ 
stat and with radioactive sources illuminating a 
12-element detector in our apparatus for neutron 
beta-decay. In addition, the multi-element de¬ 
tector allowed for an additional measurement ap¬ 
proach in which for any given crystal, a continuum 
of background gamma rays incident on the other 
11 crystals produced escaping bismuth K x-rays 
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that are detected by that given crystal. This par¬ 
ticular determination of nonproportionality at the 
Bi K x-ray energy (mean of 78.7 keV) was per¬ 
formed in a magnetic field of 4.6 T. As in prior 
published studies, we determined the nonpropor- 
tionality by comparing the center of the peak in 
the pulse height spectrum for a given gamma-ray 
energy deposit E , Ce, to that expected based on 
linear scaling from the center for the 661.7 keV 
gamma rays from a 137 Cs source, (7661.7- This def¬ 
inition applies for sufficiently high photoelectron 
yields so that a symmetric peak is obtained, which 
applies to all of our results except the 55 Fe source 
results presented at the end of the paper. For this 
definition, the measured nonproportionality Ne at 
energy E expressed in percent is given by 


N e = 100 


_ Ce _ 

(£7/661.7) 0)61.7 


% 


(i) 


Hence for a ideal scintillator, Ne is unity for 
all energies. In Sec. [2j we describe the apparatus 
we employed for our measurements of nonpropor- 
tionality. In addition to presenting our results in 
Sec. [3j we compare to previous reports. The com¬ 
bination of data taken under different conditions 
yields a consistent picture for BGO nonproportion¬ 
ality that should be useful for researchers employ¬ 
ing BGO at low gamma ray energies. In Sec. [4] we 
present results for data obtained in the test cryo¬ 
stat with 6 keV x-rays from an 55 Fe source and in 
Sec. [5] we summarize the paper. 


2. Apparatus and measurements 

2.1. Test cryostat measurements 

The test cryostat apparatus has been previously 
described in its use for studies of APD response [9]. 
For these nonproportionality studies, a 12 mm 
by 12 mm by 200 mm BGO crystal pm coupled 
to an APD with an active area of 13.5 mm by 
13.5 mm m was thermally sunk to the bottom 
of the inner reservoir of a cryostat (see Fig. [TJ. 
As discussed in Ref. [8], the opposite end of the 
BGO crystal was painted with reflective paint to 
improve the light collection efficiency. The long 
sides of the BGO crystal were covered by two lay¬ 
ers of aluminized Mylar (total of 12.4 /im Mylar 
and 0.065 /am aluminum). The long-side regions 
near the end faces («12 mm) had additional layers 
of polytetrafluoroethlyene (PTFE) tape and alu¬ 
minized Mylar for mechanical protection. (^300 
/am PTFE and 6 /im Mylar with 0.033 /im alu¬ 
minum). The cryostat’s vacuum jacket included 
a 3.5 cm diameter tube with a 7 cm diameter 
metal-sealed flange through which we illuminated 
the crystal about midway along its 200 mm length 


(see right side of Fig. [l]). The light collection ef¬ 
ficiency for gamma rays incident at different loca¬ 
tions along this crystal can vary by « 8 % from the 
center, with greater light collection efficiency ob¬ 
tained at either end of the crystal. This geometric 
effect is independent of the incident energy [8]. For 
these nonproportionality tests, the tubing for the 
vacuum port and the cryostat itself provided suf¬ 
ficient collimation to ensure that only the central 
region of the crystal was illuminated. For the high¬ 
est energy (thus most penetrating) source, this was 
confirmed by a check with additional collimation 
that consisted of 1.2 cm thick, 2.8 cm tall copper 
pieces clamped in place above and below the hor¬ 
izontal 3.5 cm tube. 

An oil-free turbomolecular pumping station was 
used to reduce the pressure in the chamber to a 
typical base pressure of on the order of 10 -6 mbar. 
Upon filling the cryostat with liquid nitrogen, the 
typical pressure was on the order of 10“ 7 mbar. 
The APD signal was sent to a charge-sensitive 
preamplifier [T2], followed by an amplifier set for a 
gain of 100 and a 6 /as shaping time constant. The 
output of the shaping amplifier was registered by a 
multichannel analyzer (MCA) [13]- The lower level 
discriminator of the MCA was typically set slightly 
below the detector noise level at ~ 10 channels, for 
a 1000-channel full range of ^1000 keV. The APD 
was operated at a typical bias of 1380 V (?^40 V 
below breakdown). Using studies with a precision 
pulser, we found that our electronics for the test 
cryostat apparatus exhibited a slight nonlinearity 
that was consistent with an offset of one channel 
out of 1000 in the MCA. To determine the BGO 
nonproportionality, we corrected for this nonlin¬ 
earity, which yielded a shift of 4.3 % for our lowest 
energy datum at 23.3 keV, decreasing to 1 % at 60 
keV. 

A thin polyimide window in the side port al¬ 
lowed the sources to be outside the cryostat vac¬ 
uum space so they could be easily manipulated. 
The gamma-ray energies and sources employed in 
the test cryostat apparatus are listed in Table [l] A 
137 Cs source was used to provide the 662 keV ref¬ 
erence. The mean energy deposit listed was deter¬ 
mined for multiple-line emissions from the follow¬ 
ing components in keV, using the weighting per¬ 
centages [IT in parentheses: 109 Cd - 22.1, 25.1 
(83.7*0.026/17.7*0.078); 133 Ba - 30.6, 31.0, 35.1, 
35.9 (54, 100, 29, 7); 153 Gd - 97.4, 103.2 (29, 21); 
and 57 Co - 122.1, 136.5 (85.5, 10.7). For the low 
energy lines of the 109 Cd source, the weighting in¬ 
cluded the factors shown to account for energy- 
dependent transmission through a 0.25 mm thick 
stainless steel cover on this particular source. To 
evaluate the nonproportionality, the pulse height 
spectrum for a given emission was fit with a sin- 
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Figure 1: Photographs of test cryostat apparatus. The left photograph shows the crystal and APD 
with the cryostat vacuum jacket and radiation shield removed. At the top of the photograph the inner 
reservoir of the dewar is visible. In the right photo (taken with cryostat rotated by 90 degrees) the 
cryostat is closed up and the vacuum port for illumination of the center of the crystal with sources 
located outside the cryostat vacuum is shown. (Color online only.) 
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E [keV] 

source 

N e [%] 

23.3 

luy Cd 

73.4 (3.0) 

31.7 

133 Ba 

81.2 (2.5) 

59.5 

241 Am 

88.8 (1.0) 

81.0 

133 Ba 

91.6 (1.5) 

88.0 

109 Cd 

90.7 (1.0) 

99.8 

153 Gd 

86.5 (1.5) 

123.7 

57 Co 

89.2 (1.5) 

185.7 

235jj 

93.3 (2.0) 

661.7 

137 Cs 

100.0 


Table 1: Mean gamma-ray energy deposits and 
nonproportionality Ne for the sources employed 
with the test cryostat apparatus. For sources with 
multiple, unresolved emission lines, the first col¬ 
umn shows the mean energy E deposit (see text). 
The numbers in parentheses show the one-sigma 
uncertainty in the last digits, which were estimated 
based primarily on reproduceability and uncer¬ 
tainty in the determination of peak locations. The 
nonproportionality is defined to be 100 % for the 
photopeak produced by the 661.7 keV gamma ray 
from 137 Cs. 

gle Gaussian, except for the multiple-line emis¬ 
sions noted above for which multiple-line Gaus- 
sians were employed. In the case of multiple-line 
emission, a single nonproportionality value was 
used for all the Gaussians and is listed for the mean 
energy deposit. 

Fig- 0 shows the background-subtracted pulse 
height distributions for three of the seven sources. 
The background was determined by counting with¬ 
out the source. However, we found that the ap¬ 
parent threshold was higher with a source than 
without a source and it varied between sources. 
The response observed near the 10 keV threshold 
is due primarily to this source-related response, 
with a smaller contribution from imperfect back¬ 
ground subtraction. 

2.2. Measurements in neutron beta-decay appara¬ 
tus 

The apparatus for neutron beta-decay has been 
described elsewhere [U 0 H]. In brief, a neutron 
beam passes through the cylindrical bore of a su¬ 
perconducting magnet and gamma rays emitted 
during neutron beta decay are detected by an an¬ 
nular array of twelve 12 mm by 12 mm by 200 mm 
BGO crystals. As in the test dewar, each crys¬ 
tal is coupled to an APD with an active area of 
13.5 mm by 13.5 mm. However, the electronic pro¬ 
cessing and spectrum analysis were different than 
that performed for the test dewar studies. The 
preamplifier output signals were directly fed into 
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Figure 2: Pulse height distributions obtained in 
the test cryostat apparatus for three of the seven 
sources employed ( 109 Cd, 133 Ba, and 137 Cs). The 
resolution is primarily determined by photoelec¬ 
tron statistics. 


two 8-channel GaGe cards [15], which were set to 
digitize the signals with 40 ns per channel resolu¬ 
tion in 80 /is long traces. The decay time constant 
of the preamplifiers is 50 /as. The amplitude of 
the digitized signals (see Ref. [8] for an example) 
were determined by fitting the signals to a tem¬ 
plate waveform. The template was a model wave¬ 
form constructed from the average of preamplifier 
signals for gamma rays with a range of energies. 

To test the detector response with radioactive 
sources, we employed a reentrant tube that allowed 
the sources to be outside the magnet’s vacuum 
chamber while illuminating the liquid nitrogen- 
cooled detector in a similar way as the photons 
from neutron decay. These tests were performed 
with the apparatus off the neutron beam line and 
with the magnetic field off. Pulse height spec¬ 
tra with 137 Cs, 133 Ba, 57 Co, and 241 Am sources 
are shown in Ref. [8], along with the results of 
modelling the detector response based, in part, 
on MCNP5 [18] predictions for the energy deposit 
spectrum in each BGO crystal. Nonproportion¬ 
ality was not included in the modeling. For this 
paper, we determined nonproportionality by com¬ 
paring the location of emission lines from sources 
to those determined from Geant4 [16] and MCNP5 
modelling with the assumption that the nonpro- 
portionality was 100 % for all energy deposits. A 
137 Cs source was used to provide the 662 keV ref¬ 
erence line for lines from a 133 Ba (nominally 31 
keV, 81 keV, and 356 keV) and a 57 Co source (122 
keV and an X-ray escape peak at 43 keV). 

An additional datum was obtained using x-ray 
escape from continuum background gamma-rays 
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E [keV] 

source 

N e [%] 

32.3 

133 Ba 

81.0 ( 2.0) 

42.7 

57 Co (escape) 

87.0 (2.0) 

78.7 

Bi escape 

94.0 (2.0) 

79.7 

133 Ba 

93.0 (1.0) 

121.9 

57 Co 

92.0 (1.0) 

356.0 

133 Ba 

98.6 (0.2) 

661.7 

137 Cs 

100.0 


Table 2: Mean gamma-ray energy deposits and 
nonproportionality Ne for the sources employed 
with the neutron beta-decay apparatus. The first 
column shows the mean energy E deposit as de¬ 
termined from our Monte Carlo modeling of the 
apparatus, without including nonproportionality. 
(Note that these values are not necessarily exactly 
the same as those in Table Q] because of a different 
experimental configuration and different methods 
for their determinations.) The nonproportionality 
is defined to be 100 % for the 661.7 keV gamma 
ray from 137 Cs. The numbers in parentheses show 
the one-sigma uncertainty in the last digits, which 
were determined from uncertainty in the determi¬ 
nation of peak locations and the variance in the 
expected peak location from the Monte Carlo for 
the 11 BGO crystals. 

detected during operation on the neutron beam 
line. Although nonproportionality may depend 
weakly on the energy deposition process [IT], in 
Ref. [4] the values obtained with escape peaks were 
the same within uncertainties as those obtained 
from direct sources. For a given BGO crystal, es¬ 
cape of Bi K x-rays from all the other crystals 
yields a peak at a mean energy of 78.7 keV. If 
the pulse height spectra from all the crystals are 
summed, this peak disappears because detection 
of the x-ray only occurs when the x-ray escaped 
from a different crystal. To improve the extrac¬ 
tion of the peak location, we combined all the data 
from the 12 crystals but also subtracted off the 
summed spectrum. In this case, the magnetic field 
was 4.6 T. The gamma-ray energies and sources 
employed in the neutron radiative decay (RDK) 
apparatus are listed in Table |2j 

3. Results 

Information on our results and those from past 
literature are summarized in Table 03 Whereas 
most results were obtained with a variety of 
radioisotopes, the work in Ref. [l] employed 
monochromatic gamma-ray beams from a syn¬ 
chrotron with resolutions varying from 2 eV at 9 
keV to 20 eV at 100 keV. In Ref. [4] a photomul- 


ID 

source 

T[ K] 

opdet 

crystal [mm] 


EL R 

295 

PMT 

not stated 


EL R 

295 

PMT 

10-20, 1-5 

1 

0,R 

295 

PMT 

9, 4 or 9 

2 

SLR 

100 

APD 

9,4 

3 

0], ER 

100 

APD 

9,4 



3, 10, 30 

PMT 

20x10x5 

4 

EL s 

295 

PMT 

5-10, 1 

5 

test, R 

90 

APD 

12x12x200 

6 

RDK, R 

90 

APD 

12x12x200 

7 

RDK, E 

90 

APD 

12x12x200 


Table 3: Identification (ID) for numerical refer¬ 
ences in Fig. [3] and information on experimental 
conditions from the literature. The information 
for our data obtained using the test cryostat (test) 
and the neutron beta-decay apparatus (RDK) are 
reported. We list the source of the gamma rays, 
the approximate temperature T of the BGO crys¬ 
tal, the method for readout of the optical photons 
(heading “opdet”), and the crystal dimensions in 
millimeters (heading “crystal”) . The source are 
synchrotron radiation (S), radioisotopes (R), es¬ 
cape peaks from radioisotopes (ER), and Bi escape 
X-rays produced from a continuum of background 
gamma rays in our neutron beta-decay apparatus 
(E). For cylindrical crystals the dimensions listed 
are in the form: diameter, length. 


tiplier tube (PMT) was used for measurements at 
room temperature and an APD PI was used for 
measurements near liquid nitrogen temperature, 
with similar results being obtained for the non¬ 
proportionality. In addition to direct gamma ray 
lines, nonproportionality was determined using es¬ 
cape peaks, again with similar results. 

Fig-i shows data for nonproportionality in the 
response of bismuth germanate (BGO). Our re¬ 
sults from the test cryostat and neutron beta- 
decay apparatus are listed in Tables [l] and [2j 
respectively. The combination of our data and 
prior published results include a variety of crystal 
sizes and sources, optical readout methods, tem¬ 
peratures, radiation sources, and analysis meth¬ 
ods. Nevertheless, all of the results generally agree 
within a few percent. The results from our test 
cryostat data are closest to the results of Ref. [1], 
whereas the results from Ref. [4, and our RDK 
apparatus are typically a few percent higher, in 
particular just below the Bi K-edge at ~90 keV. 
The single datum obtained at the bismuth K x- 
ray energy at a 4.6 T magnetic field is consistent 
within uncertainties with the nonproportionality 
values obtained at near zero magnetic field. 
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Figure 3: Data for nonproportionality in the response of BGO in the gamma-ray energy range from (a) 
0 keV to 800 keV and (b) 0 keV to 110 keV. See Tabled for identification and information on each data 
set. For clarity, the nearly continuous and precise (0.5 % uncertainties) synchrotron data from Ref. [I] 
are shown as a line with small points. The error bars for our test dewar and RDK apparatus data 
correspond to the uncertainties listed in Table [l] and Table [2j respectively. 
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Figure 4: Pulse height distribution obtained in the 
test cryostat apparatus for 1 cm 3 BGO crystal il¬ 
luminated by an 55 Fe source. The solid line is a fit 
to a gamma distribution. 


4. Nonproportionality at 6 keV 

Ref. [2] reported a monotonic increase in Ne 
with decreasing energy deposits in the range be¬ 
low the Bi L edge at ~13 keV. Whereas the ex¬ 
pected increase in the vicinity of the L-edge was 
reported in Refs. [T] and [3], Ne was found to 
continue decreasing further below the edge region. 
This discrepancy motivated us to extend our non¬ 
proportionality measurements as low as possible in 
photon energy, hence we illuminated a 1 cm 3 BGO 
crystal with 6 keV x-rays from an 55 Fe source with 
an activity of ^9 mCi. We obtained the best sig¬ 
nal to noise ratio for observing these low amplitude 
pulses by operating the APD at between 5 V and 
15 V below breakdown. For these data, the ampli¬ 
fier gain was 116, the shaping time constant was 
6 fi s, and the lower level discriminator of the MCA 
was 13 channels, for a 1000-channel full range of 
«100 keV. Under these conditions, we obtained 
the background-subtracted pulse-height spectrum 
shown in Fig. [4] To avoid extremely low channels 
on the MCA, the nonproportionality was deter¬ 
mined with respect to the 59.5 keV peak from the 
241 Am source and a correction was applied based 
on our nonproportionality value for this line. Be¬ 
cause of the peak asymmetry, we fit the data to a 
gamma distribution[20^ (see Fig. [4]). Evaluation of 
the mean value yielded 7V 6 0 = (69 =b 3) %, which 
corresponds to 1.4 times greater than that would 
be obtained based on the observed peak location. 

An asymmetric spectrum would be expected 
from Poisson statistics for a yield of only a few 
photoelectrons or less. In this case APD noise 
might mask the appearance of resolved peaks for 
specific photoelectron numbers. We modelled the 
pulse height spectrum using Poisson statistics and 
obtained the best results for a photoelectron yield 


of 2.6. However, an estimate of the photoelectron 
yield is not consistent with such a low value. For 
our large crystals we had estimated the photoelec¬ 
tron yield to be ~1000 for 662 keV gamma rays [8]. 
For 6 keV and A^.o = 100 % this estimate corre¬ 
sponds to 9 photoelectrons, but the observed non¬ 
proportionality reduces the estimated yield to 6 
photoelectrons. In addition, the small crystal em¬ 
ployed for the 55 Fe studies provided a higher yield. 
Hence the substantial asymmetry observed might 
indicate a pulse height distribution with a vari¬ 
ance greater than its mean, due to APD noise. 
The impact of the value of the yield is on the 
probability of zero photoelectrons, which yields a 
small correction to the observed nonproportional¬ 
ity. Based on our expected yield the correction is 
negligible. However, we cannot discount the possi¬ 
bility that the yield is ^2.6 photoelectrons, which 
would result in a probability of zero photoelectrons 
of exp(—2.6)=0.074. From these considerations we 
assign an asymmetric one-sigma uncertainty inter¬ 
val, hence our final value for the nonproportional¬ 
ity at 6.0 keV is 7V 6 . 0 = (69^|) %. Our result falls 
midway between the values reported in Refs, [l] 
and J3]. In Ref. [3] it was suggested that the dif¬ 
ference between their results and those of Ref. [2] 
could be related to the surface state of the inves¬ 
tigated crystals. 


5. Conclusions 

Motivated by an experiment to study neutron 
radiative beta-decay, we performed new measure¬ 
ments of the nonproportionality in the light yield 
of BGO. For the energy range from 23 keV to 662 
keV, these measurements were obtained in two dif¬ 
ferent apparatus by illuminating long crystals at 
a temperature of ~ 90 K with gamma rays from 
radioisotope sources and reading out the scintilla¬ 
tion light with APDs. In addition one datum was 
obtained in a 4.6 T magnetic field based on the bis¬ 
muth K x-ray escape peak produced by a contin¬ 
uum of background gamma rays. The combination 
of data taken under different conditions yields a 
reasonably consistent picture for BGO nonpropor¬ 
tionality. An additional measurement using 6 keV 
X-rays from an 55 Fe source was obtained with a 
small crystal. An asymmetric peak was observed, 
hence a value for the nonproportionality was ob¬ 
tained by extracting the mean value for a fit to a 
gamma distribution. 
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